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In this article we consider the three-particle decays Xbo,2 — > J/ipcc —> J/tpDD + X. We present 
the analytical formulae for the differential widths of these decays, the numerical values of their 
integrated widths and check the duality relation that connects the decays Xbo,2 '4>cc with the two- 
particle decays Xbo,2 — > J j^icc). We also study the possibility of observing the Xbo,2 — > J/ipDO + X 
mode at Tevatron and LHC colliders. 
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I. INTRODUCTION 

The bound states of heavy quarkonia (for example J /ip-, rjc- or Xf,o,2-mesons) are of great interest from both 
theoretical and experimental points of view 0. This interest is caused by the fact that these mesons can be simply 
separated at the experiment. For example, a clear signal for J/f/j-mcson (in what follows we will denote it by ip) is its 
leptonic decay ip — > e~^e~ . On the other hand the non-relativistic nature of these mesons simplifies their theoretical 
description noticeably. For the calculation of probability of the processes with such mesons one can often neglect 
the internal motion of quarks inside mesons in the hard part of the amplitude (in what follows we will refer to this 
assumption as "(5-approximation"). However the application of the (5-approximation in two-particle processes is not 
valid in some cases and can lead to large errors. The cross-section of the e~^e~ tp-qc at y^s = 10.6 GeV obtained with 
the help of this approximation for example, is about an order of magnitude smaller than the experimental value 
0. In the papers citeMa:2004qf,Bondar:2004sv,Braguta:2005gw it was shown that taking into account the internal 
motion of quarks inside mesons one increases the theoretical predictions significantly and the agreement with the 
experimental results can be achieved. The physical reason of this effect is that the intermediate particles have large 
virtuality when the ^-approximation is used, and the internal motion of quarks decreases this virtuality and hence 
raises the value of the cross section. The same enhancement can be also observed in excited charmonia production 
i. 

The method used in 0, IS can, however, lead to large uncertainties. The reason is that the quarks distribution 
functions play a significant role in this method and varying this functions we can change the result noticeably |S 0] ■ 
So, an independent way to determine these cross sections or, at least, to set some bounds on their values is desirable. 
This bounds can be determined with the help of duality relation that links the cross section of three-particle process 
e"^e~ — > tpcc with the sum of two-particle cross-sections @: 

A/th 

(icr(e+e 'pec) / + - m\ 

drUcE ^ 3 w > cr(e+e ~* M^j). (1) 

dnicc ^ 

The integration on the left hand side of this relation with respect to the invariant cc-pair mass is held over the duality 
region 

2toc < < A/th ~ 2toi3 

and the summation on the right hand side is performed over all charmonium mesons, whose masses lie in that region 
(i.e. M = ric,ip,ri'^,XcO, ■ ■ ■)■ It should be noticed that the cross-section a{e~^e~ 'pec) in the left hand side of 
the duality relation has only mild dependence on the choice of the distribution function |^ since, contrary to the 
two-particle reactions e+e~ Aiip, the virtualities of intermediate particles in this process are not fixed. 

Besides the mentioned above reactions there are also other processes where the intermediate particles have large 
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virtuality when the (5-approximation is used. In the paper [3 we considered the decays Xbo,2 ~> 'P'P and showed 
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FIG. 1: The diagrams for xj ~^ 4^cc process andx,/ g* g* vertex 

that the widths of these decays are increased when the internal motion of quarks in mesons is taken into account 
and depend strongly on the choice of the distribution functions. In the present article we will consider the reactions 
Xbo,2 ipcc that are dual to these decays and check the duality relation analogous to eq.(^ in this case. 

As it was mentioned above, when one abandon the (5-approximation, the widths of the decays Xbo,2 i^ip ^^e 
increased. In Q we have estimated the cross-section of hadronic process pp Xbo,2X ^ipX and show that it 
is quite possible to use the Xbo,2 ipip-mode for the detection of Xbj-mesons and to separate this signal from the 
background of the production of ipip pairs at Tevatron and LHC colliders. Here we will perform this analysis for the 
Xbo,2 ipcc decay. 

The rest of our paper is organized as follows. In the next section we will describe the formalism that was used to 
obtain the differential Xbo,2 "fpcc and present the analytical results. In the section 3 the duality relation is checked 
and total widths of the decays Xbo,2 — * ijjDD are given. Section 4 is devoted to the possibility of observation of these 
decays at Tevatron and LHC colliders. Finally we discuss our results. 

II- Xb0,2 i^CC 

In what follows we will use the valence approximation, i.e. restrict ourselves to the main term in the meson fock 
expansion and, in particular, neglect the contribution from the color-octet states. In this approximation it is supposed 
that the meson consists of the color-singlet quark-antiquark pair and the probabilities of such formation are determined 
from the mean values of the four-fermion operators (Oi)^ ^. A detailed description of the formalism used in our work 
can be found, for instance, in 2J. 

The diagrams of the process Xbo,2 — > V'cc in this approximation are shown on figni- The corresponding amplitudes 
can be written in the form 

M^'^ = ^{X'^X'')^^A^^iBM^uip,)ri(M^+p)l^vM), 
H1I2 ^ ^ 

where pi,2 are the momenta of final quark and antiquark, i, j are their color indices, Ui{pi) and Vj(j)2) are the 
corresponding wave functions, P and e are the momentum and polarization vector of final Vj-meson and A^J^S""^ 
denotes the amplitude of the decay of initial xtij-meson into a pair of virtual giuons with the momenta qi^2 and color 
indecies a, b (the corresponding diagrams are shown in figH^)- The widths of the Xbj-mesons decays in the considered 
mode can be expressed through these expressions with the help of well known Dalitz relation: 
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Here we use the invariant masses 



Sl,2 = M^Xi^2 = {Pl,2 + Pf, Sl2 = Mlxi2 = {pi + P2? , 

that are Hnked by the relation si + S2 + si2 = M'^ + + 2m^. The bar in eq.JSJ stands for the average over the 
polarization of the initial particle (in the case of tensor meson, J — 2) and summation over the polarization of final 
particles. With the help of this relations the following expression for the differential width of the decay Xbo V'cc 
was obtained: 

dVixba V-cc) 1310727ra4CFC 



dxidxu 9iV3(42;i - e)^iAx2 - ^2)2(4 + ^2 _ 2xi - 2x2Y 



X 



where = 3, the color factor Cp = (N^ - l)/'2Nc, S, = M^/M^, the mean values of the four-fcrmion operators 
(d)^ X can be expressed via the values of the mesons' wave functions in the origin (see next section for details) and 

the coefficients Cn^ are defined as 

C^°) = -256x1(16x1^ + 48a;i'*(-l + xiz) -96a;i2(-l + 2:12) V 
8xi3(ll- 16x12 + 5x122) + 

(5 + Xi2)^(-1 + 3X12 - 11X12^ + 9X12^) + 

xi(65 - 152x12 + 302x12^ + 72xi2^ + xi2^)), 
Cf^ = 64(288x1^ + 16xi''(-ll + 51x12) + 
32xi^(-l - 20x12 + 21x12^) + 

(5 + Xi2)^(3 + 3X12 - 23X12^ + Xi2^) + 

16x1^(30 + 201x12 + 62x12^ + 19xi2^) + 

2xi(-205 - 1296x12 + 1394x12^ + 864xi2^ + 107xi2^)), 

C^"^ = -16(-1515 + 2096x1^- 7688x12 - 3562x12^- 392x12^ + 5x12^ + 

32x1^(73 + 167x12) + 16x1^(151 + 304xi2 + 261x12^) + 

128xi(-16 + 141x12 + 92x12^ + 17xi2^)), 
C'f^ = 64(-1180 + 492x1^ -1125x12 -316x12^-11x12^ + 

2x1^(493 + 515x12) + 2xi(833 + 980xi2 + 335x12^)), 
Cf^ = -8(-3961 + 1742xi2- 3152x12 -319xi22 + xi(4258 + 2654x12)), 
Cf^ = 4(-2479 + 258x1 - 825x12), 
C("' = 1515. 

For the decay Xb2 ~^ ''pcc we have 

dT{xb2 ^cc) _ 2621447ra4CF^ 



dxidxi2 45iV3(4xi - ^2)^(4x2 - C^)2(4 + ^2 _ - 2x2^ 



{Ol)^{Ol)^ A 2n^(2) , . 



= -256xi(-l + xi+xi2)(49 + 16xi* + 32xi3(-l + xi2) + 
8xi2(l + (-20 + Xi2)xi2)- 

8xi(-l + Xi2)(l + Xi2(16 + X12)) + 

xi2(-100 + xi2(158 + Xi2(-20 + 9x12)))), 
= 64(147 + 288x1^+ 429x12 - 266x122 + 98x12^ -25xi2'' + xi2^ + 
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16a;i'^(-47 + Slxia) + 32a;i3(17 - 92x12 + 21xi2^) + 
16a;i2(-30 + 261a;i2 - 154x12^ + 19x12^) + 
a;i(694 - 3144x12 + 3460x12^ - 648x12^ + 214x12'')), 
= -16(1149 + 2096x1'' - 4148x12 + 1622x12^- 452x12^ + 5x12^' + 
32xi^(-143 + 167x12) + 16x1^(493 - 776xi2 + 261xi2^) + 
32xi(-253 + 645x12 - 172xi2^ + 68x12^)), 

CP = 64(-772 + 492x1^ + 567x12 - 184x12^- 11x12^ + 

2xi2(-383 + 515x12) + 2xi(1115 - 544xi2 + 335xi2^)), 

Cf) = -8(1937+ 1742x1^- 836x12 - 319x12^ + 2xi(-607+ 1327x12)), 

Cf^ = 4(-19 + 258x1 - 825x12), 

= 1515. 



III. NUMERICAL RESULTS AND DUALITY RELATION 



Using expressions presented above for the Xbo,2 — *■ tpcc branching fractions we can obtain their numerical values. In 
our calculations the following values for c-quark, x- and i/'-mesons were used: 

mc = 1.5 GeV, 
= 3.097 GeV, 

= (AV,+M^,J/2 = 9.886 GeV, 

and the strong coupling constant was set to be as — 0.2. The mean value of the four-fermion operator (Oi)^, in 
eqs.Q, l@J is related to the wave function of the -i/;-meson at origin or the leptonic width of this meson: 

(Oi)^, = 2N,\ibm^ = ^^^^(^ ^ e+e-)ml = 0.261 GeV^ 

Since we do not consider loop corrections we can use only leading order results in these expressions. The matrix 
element (Oi)^ can be expressed through the value of the derivative of Xb-meson wave function at the origin: 

In this work we use |i?'(0)p = 1.34 GeV^ ^3 and total widths of X6o,2-niesons in color-singlet approximation are 
equal to 

rixbo^gg) = 96a2Hl)^~ 0.544 MeV, 

In the works [lllll^ it was shown that for P-wave mesons (for example Xbj) the contribution of color-octet states is 
not suppressed with respect to singlet ones. Our estimates, however, show that this effect does not lead to noticeable 
change of Xbo,2 total widths (actually, Xbo widths is increased by approximately 5%, Xb2 width by 20%) and we will 
neglect these corrections. 

In the figI21 the distributions of Xbo ~^ "^cc (dashed line) and Xbo ~^ V^cc (solid line) versus the invariant mass of 
quark-antiquark pair TOcc are shown. For the large values of rricc the hadronization of this pair into DD is possible. 
For small values of nice the production of DD is forbidden and the quark-antiquark pair has to form a bound state, 
i.e. charmonium. This fact makes it possible to write a duality relation that is analogous to the one presented in j3|: 

A/th 

dmc,^^^^^^-l^ = Y^Brix^r ^ i^M). (5) 
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In the left hand side of this relation the integration is held over the duality interval 

2mc < Trice < -^^th — 2to£) + dM, 

and the summation on the right hand side over the (cc) mesons, whose masses belong to this interval. Because of 
charge parity conservations this should be vector mesons and excited mesons do not change the sum in the r.h.s. of 
eq.© significantly, so we will use only M = tp, ip{2S) and ■!/'(3770). We will also suppose that the excited states differ 
form ■(/'-meson only by the mass and the value of mesonic constant that can be determined from the width of the 
leptonic decay. That is why the following relation is valid: 

Br{xbj^^M) - , _/ ./ Br x&./^ V'V' ■ 

2 V{tp ^ e+e ) Mm 

According to Q the reasonable values of the parameter 5M belong to the interval 

0.5 GeV < 5M < 1 GeV, 

that is shown on figElby vertical lines. This restriction sets the bounds on the Xho.2 4''^J branching fractions: 

2.36 • 10-^ < Br (xbo ^ V'V') < 4.13 • 10-^ 
2.12 • 10""* < Br {xb2 -> ^/'V') < 3.53 • 10"". 

As in was mentioned above, the theoretical predictions for this branching fractions depend strongly on the choice 
of meson distribution functions and the presented above restrictions exclude some variants. In particular, from all 
distribution considered in only the distribution (j){x) ^ a:'^(l — xY 11^1 gives the branching fractions that are in 
these intervals 

BrCxbo^V-V-) = 2.88-10-5, Br (xfc2 ^ V'V') = 2.5 • 10"^ 
Using these numbers we obtain form eq.© the values 

5M = 0.75 GeV 

for the scalar meson XbQ and 

5M = 0.74 GeV 

for Xb2- 

At rricc > Afth the hadronization of the quark-antiquark pair into DZJ-pair and light mesons is possible. We 
will assume that the probability of this transition is 100% and that in this region the reactions Xbj ~> V'cc and 
Xbj ipDD + X are identical. The branching fractions of these decays are equal to 

Biixbo ^i^DD + X) = 1.19•10-^ 
Br {xb2 ^ ^DD + X) = 7.06-10-''. 
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IV. pp iPDD + X 

In paper we considered the possibility of observation Xb-mesons in the mode Xb V'V' Tevatron and LHC 
and separation this signal from the background of non-resonance V-'V' production. In this section we will present the 
analogous analysis for Xbj ipDD + X reaction. 

The cross section of the reaction pp ^DD + X can be expressed through the cross section of the partonic 
subprocess 



E 



dxadxbfa/p{xa)fb/p{xb)a-{ab ^pcc), (6) 







where sum is over partons a and 6, Xa.b are the longitudinal momentum fractions of these processes and fa/p and fi,/p 
are their structure functions in the initial baryons. In our calculations we used for this functions the results of the 
work |14i] . If we are interested in central region production where Xa = xi, = M^l^fs < 1 (Vs is the energy of initial 
particles in the c.m.s. frame, ^/s — 2 TeV for Tevatron and = 14 TeV for LHC), than the main contribution to 
the eq.lO gives the gluonic pair (i.e. a = b = g) and it can be written in the form 



da{pp TpDD + X) _ &gg_ 
ds s 



L. (7) 



Here s = M^XaXt = Mgg is the squared energy of the gluonic pair in c.m.s, agg ~ o^igg — > ipDD + X), and the 
dimensionless factor 

fg/p{Xa)fg/p{xb)dx (8) 

Xa+Xb 



describes the hadronic luminosity. 

The cross section of the partonic subprocess is the sum of resonance cross sections aj = cr{gg XbJ ipDD + X) 
and non-resonance background. The latter process was studied in detail in j^, where a rigorous description of 
gg — > ipcc subprocess is presented. Since we are interested only in the s dependence of the cross section of this 
process, it is convenient to use the parameterization 

that can be found in that paper. The resonance cross sections can be obtained from the presented above branching 
fractions with the help of Breit-Wigner formula: 

. ^ (2J+1) TT T{xb.j^g9)y{Xbj^^DD + X) 
16 ^Mgg-M^,,Y+TI^JA 

The widths of the resonances T^^^j are, however, small in comparison with the detector instrumental error A. We will 
take this error into account by the means of a simple substitution 

„^ ^ (2J+1) A TT V(^^,j ^ gg)T{xbj ^ ^DD + X) 

16 ^X.,M^gg {Mgg-M^,,f + ^^/A 

One can easily see this substitution does not change the value of the integrated cross section (i.e. / dMgg&j = 
J dMgg&f). In fig|3|the dependence of the total partonic cross section on the invariant mass of the gluon pair is 
shown (the value of the instrumental error A = 50 MeV was used). 

Let us now return to the cross section of the xti0,2-ineson production. Because this resonances are narrow we can 
neglect the dependence if the partonic luminosity L on s on the width of the resonance and the integral in (0 can be 
written in the simple form 

E <y{pP^XbJ+X) = _-^L^(2J+l)r(xw^ffff), 

J=0,2 X J 

where the value of the partonic luminosity L is defined according to formula JHl . For the distributions of gluons in the 
initial (anti)protons we used the parameterization presented in the work |il4 | with the value of virtuality — . 
Our calculations give L{M^) = 450 at Tevatron energies {^/s = 2 TeV) and L(M^) = 2.9 • 10^ for LHC (^s = 14 
TeV). In table I we present the cross sections of Xbo,2 productions and reactions pp XbjX t/jDD + X at these 
energies. 
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FIG. 3: The full partonic cross section (in pb) versus the gluon pair invariant mass Mgg (in GeV). The value of the instrumental 
error is A = 50 MeV. In subfigure (b) we show the region Mgg « M-^ in detail. 



J 


yfs, TeV 


a{pp Xbj + X), fih 


(^{PP Xbj +X ^ i)DD + X), nb 





2 
14 


0.25 
1.5 


0.03 
0.18 


2 


2 
14 


0.32 
2 


0.26 
1.41 



TABLE L The cross sections (t(pp xbj + X) and (t(pp — » xbj + X ^ ^pDD + X) at Tevatron and LHC 

V. CONCLUSION 

In the recent papers 0, 0, |3| it was shown that the internal motion of quarks in mesons increase the probabilities of 
two-particle reactions significantly and, in particular, an agreement between theoretical predictions and experimental 
values of the e+e~ —^ iprjc cross section and Xco ~^ 4>4' branching fraction can be reached. A poor knowledge of the 
distribution functions of quarks inside mesons can lead, however, to large uncertainties and an independent way to 
determine these probabilities or set some bounds on their values is desirable. For example, one can use the duality 
relation, that links the probabilities of three-particle reaction and two-particle processes that are dual to it. In our 
paper we have considered the three-particle process Xbo,2 — > V'cc and found the restrictions on the x&o,2 V'V' 
branching fractions: 

2.36 • 10"^ < Br (xbo ^ V'V') < 4.13 • 10-^ 
2.12 • 10"^ < Br (xb2 ^i^) < 3.53 • lO"'^. 

These restrictions, in particular, show, that from all distributions considered in only (j)(x) ~ a;^(l — a;)^ is allowed. 

In the kinematically allowed region the decays Xbo.2 ''pcc correspond to the decays x.bo,2 ipDD + X. Our 
estimates for their branching fractions are 

Bi{xbo-^i^DD + X) = 1.19•10-^ 
BT{xb2 ^ ^DD + X) ^ 7.06-10"'^. 

We have also considered the production of Xbo,2 mesons at Tevatron and LHC colliders and shown that it is possible 
to use the Xbo,2 ijjDD + X decays to detect these mesons and to separate this reaction from the background. 
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